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“Rising”  Sciences 
in  the  21st  Century 


Szczepanl990 


*  Digital  Matter 


•  Biology 


Author:  MesserWoland  and 


Medicine 


US  Air  Force 


Autonomy 


Author:  Steve  Jurvetson 
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An  Example  from  Biology 


100  Years  Ago 

-  System:  Linguistic 
Taxonomy  /  Inferred 
Heuristics 
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•  Today 

-  System: 
Hierarchical 
Informatics 


-  Tools:  Microscope, 
Cell  Culture, 
Notebooks 
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-  Tools:  Sequencers, 
Simulators,  Editors, 
“Big  Data” 
Analytics 


Which  one  does  composite  resin  chemistry  resemble? 
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The  “End  Group  Fallacy” 


*  The  “End  Group  Fallacy”  is  the  bias  that  results  from  the 

current  dominant  classification  and  inferred  heuristics  system 
for  composite  resins. 

*  Resins  are  described  according  to  a  reactive  group  class, 
“epoxy”,  “BMI”,  “polyimide”,  with  the  inference  that  key 
properties  are  shared  among  members  of  each  group.  From  the 
standpoint  of  how  to  classify  resin  systems,  such  a  system  is 
sensible. 

*  In  reality,  a  very  large  number  of  properties,  including  Tg  and 
TOS,  do  not  fall  into  separate  classes  determined  by  end  group 
chemistry.  Rather,  many  of  these  properties  are  sensitive  to  the 
topology  and  architecture  of  the  networks. 
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Structures  Architectures 


•  There  are  three  levels 


of  hierarchy  in 
networks: 


Structure  is  the  number, 
type,  and  geometric 
relationship  of  atoms  in 
repeated  groups 

Architecture  describes 
the  number,  type,  and 
geometric 
interconnection  of 
repeated  structures 

Topology  describes  the 
number,  type,  and 
geometric  relationship  of 
repeated  architectural 
units 
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Digitizing  Chemical  Structures 


Roberto  Todeschini,  Viviana  Consonni 

Handbook  of 

Molecular 

Descriptors 


Methods 
and  Principles 
in  Mcdieinal 
Chemistry 

Volume  1 1 


Edited  by 
R.  Mannhold, 
H.  Kubinyi, 

H.  Timmerman 


•  Methods  for  converting  small 

molecule  chemical  structures 
to  strings  of  information  are 
well-established  and  in  wide 
use 

•  Extensions  for  polymers  have 
been  published  but  are  not  as 
extensive 

•  Further  extensions  for 
networks  are  straightforward, 
but  not  much  used 
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*  Even  the  most  complex  cure  reactions  might  have  ~100 
repeated  chemical  structures. 

*  Architectural  features  can  be  described  in  terms  of  just  seven 
groups 

*  Topology  can  be  described  in  terms  of  a  very  few  parameters 
(typically  2  or  3) 

*  A  “cheminome”  for  a  composite  resin  is  likely  not  more  than 
10,000  units  x  10,000  unit  scale  per  each  «  106  bits 

*  In  reality,  most  properties  are  controlled  by  «  10  units. 

Modern  informatic  methods  can  determine  how  to  construct  the 
units  with  the  most  valuable  information  content. 

*  Therefore,  no  reason  that  chemometrics  /  informatics  cannot 
be  used  with  composite  resins 

AFRL: 
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Network  Automata 


Mayor  et  al.,  Macromolecules 
2011,44,  8106  (LLNL). 


•  In  correspondence  with  cellular 
automata,  a  system  of  differential 
equations  describes  the  evolution  of 
structures,  architectures,  and 
topologies  in  polymer  networks 

•  Identical  to  reaction  kinetics  at  the 
structural  level 

*  Can  include  mechanical  effects 

*  Demonstrated  for  lightly  cross-linked, 
low  Tg  networks  using  several  hundred 
thousand  units  with  sparse  topology 


Figure  3.  (a)  Representation  of  typical  bead— spring  cell  employed  in 
the  mesoscopic  model  approach.  Cross-link  junctions  are  represented 
schematically  as  beads  connected  through  springs  (straight  bars)  which 
serve  as  the  polymer  chains.  For  simplicity  of  the  image,  the  bead  radius 
has  been  chosen  arbitrarily,  (b)  Breakout  consisting  of  one  central  bead 
connected  to  nearest  and  next-nearest  neighbors.  The  dashed  line 
represents  schematically  the  real  chain  that  has  been  replaced  by  a 
single,  effective  entropic  spring  with  fc  -  (3fcsT)/((R2)). 
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•  Cure  kinetics  +  diBenedetto  equation  =  Tg  development 
kinetics  --  8  parameters  predict  Tg  of  network  through  any 
process  -  6  DSC  experiments  provide  all  needed  parameters 


Catalyst:  30  :  1  nonylphenol : 
Cu(acac)2@  2  phr  (160  ppm 
Cu) 
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V/  Example:  Bio-Based  Monomers 
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The  structure  of  bio-based  monomers 
and  networks  with  an  X-L-X 
architecture  (X  =  phenyl  cyanate  ester 
/  phenyl  cyanurate)  is  quantified  using 
1 0  parameters 
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Predictive  Models: 

Partial  Least  Squares  Approach 


#  of  Factors  #  of  Factors  ("EXP"  =  Experimental) 


Partial  least  squares  is  a  typical  informatics  technique;  it  looks  for  correlated 
parameters  that  simultaneously  describe  the  most  variation  in  both  input  and 
output  data  sets 

The  method  is  iterative;  to  generate  subsequent  regression  components,  the 
method  is  repeated  for  the  residuals  from  the  previous  set  of  predicted  and 
experimental  values 
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Predictive  Power  for 
Tg  at  Full  Cure 


"D 
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CD 
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#  of  Factors 

Error  characteristics  of  partial  least  squares  model  for  glass  transition 
temperature  at  full  cure. 

These  results  indicate  there  are  only  about  4  important  structural  parameters 
Effects  of  topology  are  not  considered 

i  5*  i 
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Example:  Bias  in  Predictive 
Parameter  Sets 
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Comparison  of  Predicted  and 
Experimental  Melting  Properties 


ASm  (kJ/mol  K,  monomer) 

BADCy 

SiMCy 

ESR255 

STT3 

ASm  (kJ/mol  K,  Yalkowsky) 

84 

84 

98 

98 

ASm(kJ/mol  K,  experiment) 

80.0  ±  1.4 

82.1  ±0.6 

75.0  ±  1.9 

74.8  ±  0.8 

ASm°  (kJ/mol  K,  Chickos) 

70 

78 

88 

95 

ASm°  (kJ/mol  K,  experiment) 

69  ±3 

81  ±1 

50  ±  14 

55  ±3 

Tm  (model  ASm  &  exp.  AHm) 

Yalkowsky  (°C) 

66 

54 

24 

24 

Chickos  (°C) 

73 

50 

42 

29 

Experiment  (°C) 

82.1  ±  0.2 

60.4  ±0.1 

115.9  ±0.2 

117.5  ±0.1 

•  Yalkowsky  model  over-predicts  entropy  of  melting  for  tricyanates,  in  part  because  the  rules 
for  counting  anisotropy  do  not  consider  star-like  arrangements,  and  a  triphenyl  substituted 
sp3  is  still  counted  as  flexible.  These  factors  explain  about  70%  of  the  error. 


'Chickos  model  has  a  similar  pattern  of  predictive  success,  perhaps  because  “bis-like” 
prolate  organic  compounds  are  more  studied  than  “tris-like”  “pitchfork”  structures 

Conversion  from  ASm°  to  ASm  based  on  Acp  m  =  const.  =  ASm ;  ASm  =  ASm°  /  [1  -  In  (Tm  /  298)  ] 
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Predictive  Parameter  Sets, 
M&S,  Databases,  and  Heuristics 


•  A  new  bit  of  data  for  a  database  is  most  valuable  when  it  describes  the 
least  known  aspect  of  the  data 

•  Good  M&S  can  derive  rules  (such  as  elastic  modulus  as  a  function  of 
topology)  for  reliably  estimating  the  lesser  known  parameters  from 
accessible  experiments 


>  Databases  and  heuristics  can  be  used  to  train  and  validate  predictive 
systems,  but  are  poor  substitutes  for  QSARs 

Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 
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Summary 


*  Significant  advances  in  materials  science  for  composites  can 
be  achieved  by  moving  from  a  “heuristic  classification” 
approach  to  an  “informatics”  approach  to  describe  chemistry 

*  Although  end-group  chemistries  play  an  important  role  as  a 
predictive  factor  in  composite  resin  processing  and 
performance,  they  exist  within  a  framework  of  structure- 
architecture-topology;  often,  other  parts  of  the  framework  are 
more  important  predictors  of  performance 


*  Even  complex  polymer  network  structures  can  be  described  in 
terms  of  a  comparatively  small  number  of  predictive  parameter 
sets;  fundamental  scientific  insight  and  access  to  a  variety  of 
chemical  structures  is  required  in  order  to  understand  what  the 
most  valuable  predictive  parameters  are 
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